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Three-dimensional Imaging System 



The current invention concerns a system for imaging simultaneously multiple layers within a 
three-dimensional object field and has applicability in fields including optical information 
storage, imaging short-timescale phenomena, microscopy, imaging three-dimensional object 
structures, passive ranging, laser beam profiling, wavefront analysis and millimetre wave 
optics. 

The use of an undistorted amplitude grating to produce identical images of a scene in several 
diffraction orders is known. Most of the energy is concentrated in the zero order with most of 
the remaining energy being contained in the +1 and -1 orders. Phase or phase and amplitude 
gratings may be used to change the distribution of energy in the different diffraction orders. 



It is also known that distortions of such a grating (i.e. dislocations in a direction perpendicular 
to the grating lines) may be used to produce phase changes in the optical system and thus shape 
the wavefront in the back focal plane of the system. This effect has been used to separate 
redundant baselines in a masked-aperture system using a dislocated grating and has formed the 
basis for computer generated holograms for many years (P M Blanchard, A H Greenaway, R N 
Anderton, R Appleby, 'Phase calibration of arrays at optical and millimetre wavelengths', J. 
Opt Soc. Am. A., Vol 13, No. 7, ppl593-1600, 1996; G Tricoles, 'Computer generated 
holograms: an historic review', Appl. Opt., Vol 26, No. 20, pp435 1-4360, 1987 and M Li, A 
Larsson, N Eriksson, M Hagberg, Continuous-level phase only computer generated holograms 
realised by dislocated binary gratings', Opt. Lett., Vol. 21, No 18, ppl516-1518, 1996). 



The imaging of a three-dimensional object using a 'through-focal series' is also known. By this 
method a sequence of images of the object are taken with the optical system focused on 
different planes in the object field. An alternative approach forms simultaneously a matrix of 
images recorded through a matrix of lenses, each of which provides a different focus condition. 
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A disadvantage of the 'through-focal series' is that because the images are recorded 
sequentially it is ill-suited to imaging the three-dimensional structure of dynamic processes. A 
disadvantage of the second approach is its complex design and that the resolution obtained is 
limited to the resolution delivered by the individual lenses in the array, the diameter of each of 
which (thus image resolution) is constrained by the space into which the array may be packed. 

The storage of data in three dimensional, optically readable, storage medium is also known (S 
Jutamulia and G M Stori, 'Three-Dimensional Optical Digital Memory', Optoelectronics - 
Devices and Technologies Vol 10, No. 3, pp343-360, 1995 and K Kobayashi and S S Kano, 
'Multi-Layered Optical Storage with Nonlinear Read/Write', Optical Review, Vol 2, No 1, 
pp20-23, 1995). These papers review the media and architecture for various three dimensional 
optical memories. 

In a high performance, near diffraction limited optical system such as a compact disk player, all 
sources of wavefront aberrations must be considered. In a standard compact disk, the data layer 
is covered with a substrate several hundred microns thick. Propagation of light through this 
substrate (which is essentially a parallel plate) introduces spherical aberration, increasing the 
spot size on the data layer and degrading resolution. This effect is overcome in conventional, 
single layer, compact disk systems by building spherical aberration correction into the objective 
lens. 

In a multi-layer optical data storage medium the degree of spherical aberration is dependent on 
the depth of the data layer in the storage medium, hence when reading from each distinct layer 
a different level of spherical aberration correction is required. An aberration corrected objective 
lens is therefore not sufficient. Several patents on multi-layer optical data storage systems, 
which rely on a moving lens to focus at different depths, have suggested ways of performing 
'active' spherical aberration correction. US 5202875 suggests using a stepped block of 
substrate material which is moved across the optical beam (using a voice coil motor) to a 
position dependent on the layer being read, such that the thickness of material that the beam 
passes through is constant. Other suggestions include a pair of prisms, one of which is 
translated, a rotating disk of variable thickness and movable compensation plates. 
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All of these approaches introduce additional moving parts and complexity into the system. 

According to this invention, an apparatus for producing simultaneously a plurality of spatially 
separated images from an object field comprises: 

an optical system arranged to produce an image associated with a first focus condition; 

a diffraction grating arranged to produce, in concert with the optical system, images associated 
with each diffraction order and 

means for detecting the images, 

wherein the optical system, diffraction grating and detecting means are located on an optical 
axis and the diffraction grating is located in a suitable grating plane and is distorted 
substantially according to a quadratic function so as to cause the images to be formed under 
various focus conditions. 

In addition, the function may be include terms to produce different amounts of spherical 
aberration in each diffraction order. This could be used to correct for spherical aberration in the 
associated optical system, for example the spherical aberration associated with layers at 
different depths in an optical data storage medium could be corrected. 

The origin of the distortion function may be displaced from the optical axis. This feature could 
be used to cause alignment along the optic axis of the images associated with each diffraction 
order. 

The grating could be a single grating or could be constructed from a combination of two or 
more gratings. 

A variety of grating types could be used in the current invention for example, amplitude only; 
phase only; a combination of amplitude and phase; polarisation sensitive; a combination of two 
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gratings sensitive to different polarisations; programmable; reflective; transmissive; two-level 
(binary); multi-level (digitised) or continuous-level (analogue). 

By one implementation, images of loci of different depth within the object field (i.e. of 
different object planes) can be produced simultaneously on a single image plane. This would 
allow simultaneous reading of optical data stored in different layers within a suitable storage 
medium (e.g a multi-layer CD-ROM). 

By an alternative implementation, images of a single plane in the object field (i.e. of a single 
object plane) can be produced simultaneously on a number of different image planes. This 
would allow focusing of a single illumination source on a number of different planes. 

The above two implementations can be combined to provide a system in which the object 
planes of the former are coincident with the image planes of the latter and hence are suitably 
illuminated to facilitate imaging. 

In a preferred embodiment the object planes contain data storage elements and can be imaged 
and, more preferentially, can be illuminated by the current invention.This embodiment could be 
used to read data from a three dimensional optical storage medium. 

Another embodiment includes a dispersive system for introducing an offset to an input beam of 
radiation without altering the direction of the beam, said offset being perpendicular to the 
optical axis and proportional to the wavelength of the input radiation. This causes the 
diffraction angles at all angles to be substantially equal. 

The invention utilises a single lens or multiple lens system with a distorted diffraction grating 
to produce simultaneously a set of images of the object field in which each image in the set can 
correspond to an image of the object field recorded under different focus conditions but in 
which the full diameter of the lens system is exploited in each image in the set. For each image 
in the set, the resolution, magnification and depth of focus is that which would have been 
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obtained if a through focal series had been produced by varying only the focal length of the lens 
system. 

The grating used can be a single distorted diffraction grating or a series of such gratings. The 
gratings used may be produced by computer-generated (digitised in space and/or in amplitude) 
or by analogue (e.g. interferometric) means. 

In the following descriptions detector means a detection means comprising a spatially-resolving 
system such as a pixellated array of detector elements e.g. a charge coupled device (CCD). For 
applications where detection of the presence or absence of unresolved targets is required, the 
detector may comprise suitably-positioned, isolated detector elements. 

The invention will now be described with reference to the following figures in which: 

figure 1 shows schematically suitable grating planes, by way of illustration only. 

figure 2(a) shows schematically a conventional, undistorted, amplitude-only diffraction grating 
used in an imaging system and figure 2(b) shows the normalised intensity cross-section of the 
zero, +1 and -1 diffraction order images of a point source produced when such a grating is 
inserted in a suitable grating plane of an imaging system; 

figures 3a and 3b illustrate respectively a section of an undistorted grating and the distortion of 
a grating by a fixed amount, A; 

figure 4 shows schematically a simple imaging system of the current invention; 

figure 5 shows two computer generated amplitude gratings with quadratic distortion; 

figure 6 shows computer simulations of the intensity cross-sections of images of a point source, 
formed in the +1 , 0 and -1 diffraction orders by the current invention; 
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figure 7 shows experimental measurements of the images of a resolution target formed in the 
+1,0 and -1 diffraction orders by the current invention; 

figure 8(a) shows crossed amplitude gratings, figure 8(b) shows the defocus states of the 
corresponding diffraction orders and figure 8(c) shows a computer simulation of the images of 
a point source through the grating structure of figure 8(a); 

figure 9 shows experimental images of a resolution target obtained at different detector 
positions using crossed amplitude gratings; 

figure 10 shows schematically a simple imaging system of the current invention used to 
produce in-focus images of different object planes at a single detector plane; 

figure 1 1 shows images of an object plane containing a slide of the letter 'B' obtained by 
experiment using a quadratically distorted amplitude grating; 

figure 12 shows schematically the apparatus used to record the experimental data shown in 
figure 13; 

figure 13 shows simultaneous side-by-side images of 3 objects located in different object 
planes; 

figure 14 shows how the invention may be adapted for reading of data stored in a three 
dimensional optical storage medium; 

figure 15 shows a comparison of phase profiles and grating structures for gratings with defocus, 
spherical aberration and both built in; 



figure 16 shows qualitatively image cross-sections associated with -1,0 and +1 diffraction 
orders through a grating with 0 C 40 =1A,; 
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figure 17 shows intensity cross-sections of images of point sources on layers corresponding to 
7, 6 and 5 of figure 14 with defocus only (left hand column) and defocus and spherical 
aberration (right hand column) corrected gratings; 

figure 1 8 shows a quadratically distorted diffraction grating with 0 C 20 = 4 A. and an offset of the 
quadratic function relative to the optical axis of x 0 =^R 2 /(2 d 0 C 20 ); 

figure 19 shows an optical system for obtaining images of 2 or more object planes with equal 
magnifications; 

figure 20 shows a quadratically distorted diffraction grating and the beam positions at 3 discrete 
wavelengths needed to equalise the angles of diffraction; 

figure 21 is a schematic diagram of a dispersive system for introducing wavelength-dependent 
beam shift prior to the distorted grating; 

figure 22(a) illustrates the desired illumination of multiple layers in an object field and figure 
22(b) illustrates illumination of such layers using gratings of the type shown in figure 2 and 

figure 23 is a schematic diagram of a multi-layer optical data storage read head; 

Although the following examples relate to application of the invention in the field of optics, 
this should not be seen as limiting as the general principles of the invention are applicable to 
other wavelengths of electromagnetic radiation. 

Grating Location 

In a system in which converging and/or diverging beams exist, a suitable grating plane would 
be any plane that is normal to the optical axis and close to a lens other than a lens positioned in 
the image or object field, for example plane PI in figure la. In a system where a collimated 
beam is produced, a suitable grating plane would be any plane that is normal to the optical axis 
of the system and in the region in which the beam is collimated, for example anywhere between 
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planes PI and P2 in figure lb, or a plane described as a suitable grating plane for a system with 
converging or diverging beams. 



Grating Design 

The design of distorted gratings which might typically be used in the current invention will 
first be described. 

A standard diffraction grating consists of alternate regularly spaced strips of different 
transmissivity, reflectivity or optical thickness. When the grating is used within an imaging 
system, multiple diffraction orders appear in the image plane in addition to the unscattered zero 
order. Each diffraction order contains the same information about the object field as the zero 
order, though with different levels of intensity dependent on details of the grating construction. 
Figure 2 shows, as an example, an amplitude grating and the images of a point object formed in 
the -1, 0 and + 1 diffraction orders, (each order normalised to unit intensity) 

If the grating geometry is distorted locally, by a displacement of the strips in a direction 
perpendicular to their long axis, a phase shift is introduced in the wavefront scattered from the 
distorted region, the level of which is dependent on the amount of local distortion of the grating 
relative to its undistorted form. The level of local phase shift is related to the distortion of the 
grating through equation 1, 

, 27imA ^ 
<P = — - — Equation 1 

d 

where d is the grating period, m is the diffraction order into which the wavefront is scattered 
and A is the distortion of the grating strips relative to their undistorted position, as shown in 
figure 3. Such a distortion of the grating produces phase shifts of equal magnitude but opposite 
sign in the wavefronts scattered into +1 and -1 diffraction orders and leaves the unscattered 
wavefront in the zero order unaltered. 



It is important to note that this technique allows continuous phase values to be encoded using a 
binary (two level) grating, although the invention can also be applied to multiple or continuous 
level gratings. 
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For applications using computer-generated holograms, the distorted grating can be designed by 
dividing the grating area into a number of cells, which can be of any space-filling shape, and 
calculating the degree of distortion to be applied to the grating for each cell individually. 
Alternatively the distortion can be applied to the grating as a whole leading to continuously 
distorted strips. Both of these approaches can be implemented using computer design followed 
by grating fabrication or by using an electrically addressed liquid crystal or other electro-optic 
device. 

For non-digital production methods an alternative technique is to record holographically the 
distorted fringe pattern into an optically sensitive medium, or to use an optically programmable 
liquid crystal device to allow the grating to be changed in real-time. 

The above descriptions refer to arbitrary distortions that could be used to generate arbitrary 
phase changes on the wavefront scattered into a selected diffraction order. 

Below are described the grating distortions required to produce the defocus effects required for 
implementation of this invention. 

Defocus Gratings 

A defocused optical system has a phase shift which, compared to an in-focus image, can be 
represented by a quadratic function of the distance from the optical axis and measured relative 
to the Gaussian reference sphere (e.g. section 5.1, Principles of Optics, Born & Wolf, 
Pergammon, Edition 6, Oxford, 1980). This invention relates to a diffraction grating distorted 
as a quadratic function of distance from the optical axis of the system according to, 

A(x,y) = -^^(x 2 + y 2 ) Equation 2 

XR 

where A(x,y) is a distortion in a direction perpendicular to the grating lines (figure 3), x and y 
are Cartesian co-ordinates relative to an origin on the optical axis in the plane of the grating, d 
is the grating period, X is the optical wavelength, 0 C 20 is the degree of defocus introduced into 
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the image formed in the +1 diffraction order ( 0 C 20 ^0) and R is the radius of the grating aperture 
which is centred on the optical axis. In equation 2 a circular aperture has been assumed, but the 
invention can be applied to an aperture of any shape. 0 C 20 is the wavefront coefficient of 
defocus of the grating (the traditional defocus aberration constant equivalent to the pathlength 
difference introduced at the edge of the aperture between, in this case, the wavefront scattered 
into the +1 diffraction order and the Gaussian reference surface for that diffraction order {e.g. 
section 15-5, Geometrical and Physical Optics, R S Longhurst, Longman, Edition 3, London, 
1973}). The phase change imposed on the wavefronts scattered into the various diffraction 
orders can be calculated by combining equation 1 and equation 2 to give, 



(x,y) = m ^£* (x * + 3) Equation 3 

XR 



The quadratic phase function (equation 3) imparts a phase delay on wavefronts scattered into 
the non-zero diffraction orders such that the wavefront curvature is altered. The grating 
therefore has focussing power in the non-zero orders, and an equivalent focal length (f m ) can be 
defined for these orders, 



f - r2 -(" <Ao) 2 



Equation 4 
2m 0 C 20 



In practice it is more useful to use such a grating as part of an optical system which provides 
the majority focussing power with the grating effectively modifying the focal length of the 
lens in each diffraction order. For example, when a quadratically distorted grating is placed in 
contact with a lens of focal length f, the focal length of the combination in each diffraction 
order is given by, 



4- f (R 2 -m 2 qC^o) _ . . 

f m = — — % 2~ Equation 5 

R + 2f m 0 C 20 -m 0 C 20 



which can be approximated by 
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Equation 6 



R +2f m 0 C 20 

when R 2 »m 2 0 C 20 . The magnitude and sign of the defocus is dependent on the diffraction 
order (m). Thus a series of images of the object field with differing defocus conditions is 
produced simultaneously and side-by-side on the detector in the different diffraction orders. 

The principle of the invention can be demonstrated with reference to the -1 , 0 and +1 
diffraction orders. Referring to figure 4, the defocusing effect of a quadratically distorted 
grating can be demonstrated using an optical system (1), designed and arranged to image an 
object (2) on the optical axis (3) onto detector plane B at the normal focal plane of the optical 



A quadratically distorted diffraction grating (4) which is added to the optical system (1) 
produces two additional images of the object (2) in plane B in its +1 and -1 diffraction orders. 
In the normal focal plane B the zero order image remains in focus, whilst the images in the +1 
and -1 diffraction orders undergo defocus of equal magnitude but opposite sign. If the detector 
is moved along the optical axis either side of plane B, a plane can be reached where the 
physical defocus cancels out the defocus introduced by the grating into the diffraction orders. In 
this way the images in the +1 and -1 diffraction orders can be brought in to focus (planes A and 



The separation 8j of the image planes A, B and C is determined by the grating distortion, the 
radius of the grating aperture and the optical system through, 

^ m - 2y2m ° C2 ° - Equation 7 

2vm 0 C 20 + R 

where R is the grating aperture radius, v is the distance from the normal image plane (B) to the 
secondary principle plane of the optical system, and the approximation R»m 0 C 2 o (R is much 
greater than m 0 C 20 ) has been made. Note that if a grating is designed with defocus represented 
by (Ao = then the +1 diffraction order undergoes a defocus equivalent to nX, the -1 



system. 



C). 
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diffraction order will undergo a defocus equivalent to -nk and, through equation 7, planes A 
and C will be located either side of and at different distances from plane B. 

In the case where 2vm 0 C 20 « R 2 , equation 7 can be approximated by, 



8; «-2m — 



0 C 20 Equation 8 



and planes A and C are symmetrically placed about plane B. 

Equation 7 can be rearranged in terms of the grating defocus ( 0 C 20 ) needed to generate the 
required image plane separation (8 { ) between in-focus images in the zero and +1 diffraction 
orders (m=l), 



° C20= 2v( R v! i 5 i ) Equation 9 



Figures 5(a) and 5(b) each show examples of gratings with spherical apertures and R=12.5d, 
distorted as a quadratic function of distance from the centre, to give different levels of defocus, 
for figure 5(a), 0 C 2 o=k and for figure 5(b), 0 C 20 =2A,. These represent two examples of many 
possible grating structures and were designed by computer as binary amplitude gratings using a 
square design cell. The distinct steps in the black grating line edges are due to the X/20 
accuracy to which these gratings have been plotted. The grating lines associated with 
quadratically distorted gratings are arcs of concentric circles of different thickness and could be 
plotted as such. 

The defocusing property of these gratings has been verified with computer simulations and 
experimentally as will be described later. 

Imaging a single object through a quadratically distorted grating - Computer Simulations 
Computer simulations were performed using software written in Fortran using a Fast Fourier 
Transform (FFT) routine [Subroutine fourn from 'Numerical recipes in Fortran', W H Press, S 
A Teukolsky, W T Vetterling, B P Flannery, Cambridge University Press, 1992]. Images were 
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calculated by multiplying the FFT of the object by the optical transfer function of the grating, 
followed by an inverse FFT to generate the image. The optical transfer function of the grating 
was calculated from the autocorrelation of the grating, obtained via a double FFT technique 
using the Wiener-Khintchine Theorem [Fourier Optics : An Introduction, E G Steward, 2 nd 
edition, p95, J Wiley & Sons.]. This approach represents a simulation of incoherent imaging. 
Figure 6 shows intensity cross-sections of simulated images of a point source through a 
distorted amplitude grating designed with 0 C 20 = X (figure 5a), with the detector placed at 
planes A, B and C (figure 4). Using an amplitude grating, the intensities of the first order 
diffraction spots would, in practice, be lower than that of the zero order. In this figure the power 
in in-focus diffraction orders has been normalised to unity to aid observation. With the detector 
in plane B, the zero order is in focus and the +1 and -1 diffraction orders have defocuses of 
+1A, and -1 A, respectively. By moving the detector either side of this plane, the +1 and -1 
diffraction orders can be brought into focus. This demonstrates that the mask is generating a 
true defocus. 

Imaging a single object through a quadratically distorted grating - Experimental Results 
In order to verify the computer simulations, a grating was fabricated by photographically 
reducing an enlarged black and white picture of the appropriate pattern on to a 35mm slide. 
This provided a grating with a circular aperture of diameter 1cm, 0 C 2 o == ^ and a grating period of 
400 x 10 " 6 m (400^m). The optical system comprised two lenses with focal lengths of 50cm and 
100cm, separated by 5cm. The object, a standard resolution target, was placed one focal length 
(50cm) in front of the first lens and the detector was placed one focal length (100cm) behind 
the second lens. A white light source was used to illuminate the object in transmission and the 
grating was placed between the two lenses in the region where the light was collimated. A filter 
with a bandpass of lOnm, centred at 650nm, was placed in front of the CCD detector used to 
record the image. 

These parameters lead to an axial focal shift of -4.9cm and +5. 5cm in the +1 and -1 diffraction 
orders respectively (equation 7). Figure 7 shows the images obtained upon location of the 
detector at positions corresponding to planes A, B and C of figure 4. The figure shows the raw 
images captured by the detector and the same images after processing to increase the intensities 
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of the +1 and -1 diffraction orders (normalised), to aid observation. It can be seen that the -1, 0 
and +1 diffraction orders are brought into focus as the detector is scanned along the optical axis 
of the system. At these positions, the physical defocus is cancelling out the wavefront 
deformations introduced by the gratings, that is the grating is introducing the quadratic 
variation of phase (defocus) predicted. 

Imaging a single object through two crossed quadratically distorted gratings 
The techniques described so far can be extended by using multiple crossed gratings. If two 
gratings are crossed at right angles, the central nine diffraction orders can be usefully used. If 
the defocuses ( 0 C 20 ) of the two crossed grating are chosen to be aX and bX then, for I a-b\ 
b, the nine images of the scene that are formed in parallel correspond to nine different defocus 
conditions. Figure 8a shows an example of two crossed gratings having defocuses of 0 C 2 o == 0.5X 
and 0 C 2 o=1.5A., figure 8b shows the relative defocuses of the central nine diffraction orders and 
figure 8c shows a computer simulation of the image of a point source through the gratings 
(normalised). The image of the object in each diffraction order can be brought separately into 
focus by movement of the detector along the axis. 

The crossed grating technique has been tested experimentally using the arrangement previously 
described and the crossed grating structure illustrated in figure 8a. The experimental results in 
figure 9 show a selection of the nine images, brought into focus by moving the detector along 
the optical axis. The full range of motion of the detector between the two extremes in figure 9 
was approximately 20cm. 



Imaging multiple object planes through a quadratically distorted grating 
The function of the defocus grating can be considered in a different way. Referring to figure 10, 
if the detector is placed at image plane B, then the three images formed correspond to in-focus 
images of three different object planes 5, 6 and 7. The zero order will be the sum of the out-of- 
focus images of objects 5 and 7 and an in-focus image of object 6. If the degree of defocus is 
sufficient, a good image of object 6 will result. Similarly, objects 5 and 7 are discernible in the 
+1 and -1 diffraction orders. The grating therefore generates, side-by-side, simultaneous images 
of three different object planes at a single detector plane. The separation (8 0 ) of the object 
planes imaged in plane B is determined by the grating distortion, the radius of the grating 
aperture and the optical system through, 
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Equation 10 



2um 0 C 20 +R 2 



where 0 C 20 is the wavefront coefficient of defocus of the grating for the +1 diffraction order, R 
is the grating aperture radius, m is the diffraction order, u is the distance from the central object 
plane to the primary principle plane of the optical system, and the approximation R»m 0 C 20 
has been made. 

The resolution in depth, in terms of the minimum separation of planes in the object field that 
can be individually imaged, is dependent on the depth of focus of the optical system being 
used. The image quality obtained when using a distorted diffraction grating to image multiple 
planes within the object field will be the same as if a 'through focal series' were obtained by 
adjusting the optical system to adjust its focus to image the same planes. 

The fact that different object planes are imaged into different diffraction orders was first 
observed using a fixed detector and a single moveable object. Figure 1 1 shows images obtained 
on locating a slide of the letter 'ET of height 1mm (the object) at positions corresponding to 
object planes 5, 6 and 7 of figure 10. The optical system comprised a single lens with a focal 
length of 12cm, a CCD detector fixed at 32cm from the lens (plane B) and a grating with an 
aperture diameter of 2cm, 0 C 20 =10A. and a period of 100 x 10" 6 m (1 OO^im). Object plane 6 
corresponded to a plane 19cm from the lens and the grating parameters produced in-focus 
images in its -1 and +1 diffraction orders of object planes displaced by approximately ±4.8mm 
relative to plane 6. 

In order to demonstrate the simultaneous imaging of three object planes, the same optical 
system was used with three different objects. Three slides containing the letters 4 A' 'B' and f C 
were placed at planes 5, 6 and 7 respectively as shown schematically in figure 12. Plane 6 
corresponded to a plane 19cm from the lens LI and planes 5 and 7 corresponded to object 
planes displaced by +4.7mm and -4.9mm relative to plane 6. A white light source (not shown) 
was used to illuminate the objects in transmission, and a filter 9 centred at a wavelength of 
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650nm, with a bandpass of lOnm was located over the detector aperture. The detector (8) was 
also moved along the optical axis to the positions where the diffraction orders imaged different 
objects, corresponding to planes A and C in figure 4. These positions were approximately 
±14mm from the image plane B. 

Figure 13 shows the images recorded on the detector at planes A, B and C of figure 12. At 
position B, the zero order images object plane 6 (the letter 'B'), the +1 diffraction order 
images object plane 7 (the letter 'C') and the -1 diffraction order images object plane 5 (the 
letter 'A') of figure 12. These images demonstrate experimentally that three object planes can 
be imaged simultaneously and side-by-side on a single detector. With two crossed gratings it is 
possible to image simultaneously nine object planes side-by-side on a single detector. 

Phase Gratings 

If an amplitude grating is used as a defocusing element, the zero order is always brighter than 
the +1, -1 and higher diffraction orders. The distribution of energy can be adjusted using a 
phase grating with two phase levels. For example, a phase step of n radians can completely 
eliminate the zero order, whilst putting more power into the +1 and -1 orders, or the phase step 
can be adjusted to 0.63971 to place equal power into the 0, +1 and -1 orders. With crossed 
gratings the use of more phase steps or combined amplitude and phase gratings can be used to 
adjust the intensities of the multiple diffraction orders. 

Polarisation-Sensitive Gratings 

A distorted grating that is fabricated such that it operates only on one polarisation state of the 
incident wavefront will produce polarised, defocused images of an object in the +1 and -1 
diffraction orders and an unpolarised image in the zero order. If a second grating, that operates 
only on the orthogonal polarisation, is crossed with the first grating then two further diffraction 
orders will be produced, polarised in the orthogonal sense and displaced from the set produced 
by the first grating. Because the gratings are polarisation sensitive there is no crosstalk between 
the gratings and thus no diffraction orders are produced other than would be produced by each 
grating acting alone. If the system is required for polarimetric studies the defocus can be chosen 
to be the same for each polarisation state. Because the images are produced simultaneously the 
system is suitable for polarimetric studies of dynamically-changing scenes. 
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Three Dimensional Optical Data Storage, 

Referring to figure 14, apparatus of the invention, adapted for reading data stored in a three 
dimensional optical storage medium 1 1 is generally designated 10. The storage medium 1 1 
comprises discrete optically readable planes 5, 6, 7 having individual data storage elements (not ' 
shown) located thereon and is illuminated by means not shown. The elements could, for 
example, be an area of the plane in which a hole or pit is used to designate a logical 1 and the 
absence of a hole or pit is used to designate a logical 0. However this example is not limiting: 
other embodiments are known to, or would be obvious to, persons skilled in the art including 
multi-level systems in which the elements may exist in more than the two states of a binary 
system. The data storage elements are imaged simultaneously at detectors 12. Detectors 12 are 
capable of producing a signal dependent on the state of the storage element and could be a 
photodiode or a photo transistor. 

In order to facilitate interrogation of different data elements within each plane, the apparatus 
includes means (not shown) for effecting relative movement, in a direction perpendicular to 
optical axis 3, between the storage medium 8 and the rest of the apparatus. Such means might 
comprise an electromechanical arrangement known to a person skilled in the art. 

Inclusion of Spherical Aberration Distortion Function 

The phase profile of a dififractive element that imparts a spherical phase delay can be written as, 



where f is the element focal length (or radius of curvature introduced) and r is a radial co- 
ordinate. This expression can be expanded as , 




Equation 1 1 



2k r r 




Equation 12 
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By using diffraction gratings that are quadratically distorted only the first term of this equation 
is being considered. The omission of higher order terms is equivalent to the standard paraxial 
approximation in which the inclination of a ray to the optical axis is assumed to be small ( 
K<f). In situations where this approximation is not valid it is appropriate to include in our 
grating distortion function terms in higher powers of r 2 , which modify the quadratic grating 
distortion. In general the grating will remain substantially quadratically distorted, with the 
higher order terms acting as corrections to the grating distortion. The higher order terms in 
equation 12 are referred to here as spherical aberration terms and will be included in the 
wavefront aberration function (W SA ), written as, 



W SA (r)- 0 C 40 r 4 + 0 C 60 r 6 + 0 C 80 r* + ... = X 0 C 2j . 0 r 2j Equation 13 

primary secondary 



where r is the distance from the optical axis. 

The ability to generate different levels of spherical aberration in each diffraction order can be 
used not only to adjust the grating structure in the non-paraxial case, but also to correct for 
spherical aberration introduced by the object field, for example in multi-layer optical data 
storage systems. 

For a beam entering a parallel plate of refractive index n, with zero tilt, the primary and 
secondary coefficients of spherical aberration are given by, 

0C40 = ( " 7 1) t < NA > 4 oCeo = ^-^t(NA) 6 Equation 14 

8n 16n 

where t is the depth at which the beam is focused and NA is the numerical aperture of the beam 
(J Braat, 'Analytical expressions for the wave-front aberration coefficients of a tilted plan- 
parallel plate', Applied Optics. Vol.36, No.32, 8459,1998.). The linear dependence of spherical 
aberration on depth (t) means that all of the terms in equation 14 can be corrected with a 
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suitably distorted diffraction grating. In the following discussion only the 0 C 40 r 4 term of 
equation 14 is considered, for clarity. 



To eliminate spherical aberration the optical system 1 (without grating) of figure 14 can be 
designed to correct for the spherical aberration introduced on the beam focused on layer 6, as in 
a standard 1 -layer compact disk read head. A diffraction grating distorted according to r 4 can 
then be designed with a wavefront coefficient of spherical aberration (oC 40 ), such that the -1 
diffraction order (associated with a spherical aberration correction of- 0 C 40 ) is spherical 
aberration corrected for layer 7, the +1 diffraction order (associated with a spherical aberration 
correction of + 0 C 40 ) is corrected for layer 5 and the zero order remains spherical aberration 
corrected for layer 6. 



Table 1 shows the defocus and spherical aberration (first order only) terms, associated with 
each layer 5, 6 and 7 introduced into grating 4. 



Table 1 



Layer 


Diffraction 


Grating Correction 


Grating Correction 




Order 


(Defocus) 


(Spherical Aberration) 


5 


-1 


"(Ao 


-<Ao 


6 


0 


0 


0 


7 


+1 


^0^20 


+ oC 4 o 



The complete grating must incorporate both defocus and spherical aberration correction and 
will be distorted according to, 



A(x,y) = ^ 



(^x o ) J + (y-y 0 ) 2 (x 2 +y 2 ) 2 , (x 2 +y 2 ) 3 h5otar4mW _ 
0^20 + — 0^40 + ~ 0^60 + higher order terms 



R 4 



R 6 



Equation 15 
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where D(x,y) is a distortion in a direction perpendicular to the grating lines, d is the grating 
period, R is the radius of the grating aperture and x and y are Cartesian coordinates relative to 
an origin on the optical axis in the plane of the grating. The grating distortion D(x,y) introduces 
a phase shift onto the wavefront scattered from the grating into the m th diffraction order 
according to equation 1. x 0 and y 0 represent an offset of the origin of the quadratic distortion 
function from the optical axis of the system. Exploitation of such an offset is described later. 
Offsets of the spherical aberration terms in equation 1 5 ( 0 C 40 , 0Q0 ^ higher terms) can not be 
exploited in the same manner and are not included in equation 15. Such distortion functions 
should be centred on the optical axis. 



As an example, the level of spherical aberration to be corrected can be estimated by assuming a 
multi-layer optical storage medium with a refractive index of 1 .5806 at a wavelength of 650nm, 
a numerical aperture of 0.60 and layer separations of lOOjam. These parameters give additional 
spherical aberration of 0 C 40 =0.95A. and 0 C 60 =0.24X from layer to layer (equation 15). (J Braat, 
'Analytical expressions for the wave-front aberration coefficients of a tilted plan-parallel plate', 
Applied Optics. Vol.36, No.32, 8459,1998.). 



Figure 1 5 compares the phase profiles and grating structures for gratings with defocus only (a), 
spherical aberration only (b) and defocus and spherical aberration (c) built in. The grating lines 
in (a) are arcs of circles, whereas those in (b) and (c) are not. 

The ability of a grating distorted according to r 4 to generate spherical aberration of equal 
magnitude but opposite sign in the +1 and -1 diffraction orders has been demonstrated in 
computer simulations. The simulation modelled a wavefront with specified levels of spherical 
aberration, incident on grating (b) of Figure 15. Figure 16 shows that when a wavefront 
aberration of one wave of spherical aberration is used, the image in the +1 diffraction order is 
corrected (the image is a diffraction limited spot), whereas when a wavefront with one wave of 
spherical aberration of opposite sign is used, the image in the -1 diffraction order is corrected. 



The Y (vertical) axes in figures 16 and 17 represent intensity and the X (horizontal) axes 
represent distance along a line through the image formed. 
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The operation of a grating with built-in defocus and spherical aberration correction has also 
been demonstrated with computer simulations. With reference to Figure 14, consider an optical 
system focused and spherical aberration corrected for layer 6. Layers 7 and 5 are located 
approximately 100 x 1 O^m either side of layer 6 r associated with additional defocuses of, say, 
plus and minus one wave and additional spherical aberration of, say, plus and minus one wave. 
Point sources are located on each of the three planes. A grating with a purely quadratic 
distortion function ( 0 C 2 o=lA,, (Acr^X designed to focus the diffraction orders on layers 5, 6 and 
7 (in the +1 , 0 and -1 diffraction orders respectively) produces the image cross-section shown 
in the left hand column of Figure 17. The images of the point sources on layers 7 and 5 are 
enlarged and reduced in intensity due to the uncorrected spherical aberration. A grating 
designed with defocus and spherical aberration correction ( 0 C 2 o~lK 0 C 40 =1X) produces the. 
image cross-section shown in the right hand column of Figure 17. The spherical aberration is 
now corrected for all data layers and the images of the point sources are in focus and of 
diffraction limited size. 

Shift of Quadratic Distortion Origin 

Consider an undistorted grating ( 0 C 20 - (Ao = <Ao- °) consisting of parallel strips of different 
transmissivity, reflectivity or optical thickness. The y-axis is defined to be parallel to the strips 
in the grating and the x-axis to be perpendicular to the strips. A plane wavefront incident 
normally on the grating is diffracted into orders at angles G m to the optical axis where, 

mX 

sin8 m =— Equation 16 

d 

The angular deflection of each order is equivalent to introducing onto the incident wavefront a 
phase tilt across the grating plane of, 



<t>(x,y) = — - — . 
d 



Equation 1 7 
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When 0 C 20 is non-zero the expansion of the defocus (first) term in equation 1 5 produces a phase 
shift given by, 



<Kx,y) = 0 C 20 



27lm (x 2 i y 2 ) - 4nmx o x _ 47tm yo , 2nm(x 0 + y 0 ) 
R 2 * R 2 * R 2 X R 2 X 



Equation 18 



The quadratic dependence of equation 1 8 (all within the first term) is identical to equation 3 
(the defocus phase term obtained from a grating distorted according to a quadratic function 
centred on the optical axis). Shifting the origin of the quadratic function therefore has no effect 
on the level of defocus introduced into each diffraction order. 

The final term in equation 18 represents a dc phase offset which has no effect on the wavefront 
shape diffracted into each order. 

The second term in equation 1 8, 

<Kx,y) = 0 C 20 x 0 x Equation 19 

A»R 

represents a linear increase in phase across the x-axis of the grating plane. This phase tilt has 
the effect of changing the separation of the grating diffraction orders, whilst leaving the 
position of the zero order (m=0) and level of defocus unchanged. If x 0 is chosen to have a value 
of, 

x o = ^ . r> — Equation 20 

2d 0 C 20 

then equation 19 becomes equal to, 



i, v —2m 

<l>(x,y) = — - 



Equation 21 
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This phase tilt is equal and opposite in sign to that introduced into light scattered from the 
gating structure. The cancellation of these two terms (equation 17 and equation 21) for each 
order causes the diffraction order spacing to become zero in a plane perpendicular to the optical 
axis, corresponding to all of the orders lying along the optical axis of the system. The 
diffraction orders remain spatially separated along the optical axis according to the level of 
defocus ( 0 C 20 ) within the distorted grating. The shift of the origin of the quadratic function to 
that defined in equation 20 leads to a grating structure consisting of circular fringes, as shown 
in Figure 18. This form of grating is used as part of an illumination system described below. 

The third term in equation 1 8, 

4>(x,y) = — ^ 0 C 20 y 0 y Equation 22 

XK 

represents a linear increase in phase across the y-axis of the grating plane. This causes the 
diffraction orders (other than m=0) to move along the y-direction in the image plane. 

Through choice of Xq and y 0 , the position of a particular diffraction order in the image plane can 
be controlled, whilst leaving the level of defocus and spherical aberration unchanged. 

System with Equal Image Magnifications 

Consider a system of the type shown in Figure 10 in which multiple object planes are imaged 
onto a single image plane. In general, according to equation 10, the object planes imaged into 
the +1 and -1 orders (7 and 5 respectively) are not symmetrically located about the plane (6) 
imaged into the zero order. The magnification associated with each image (M m ) calculated 
simply as the ratio of the image distance (v) to object distance is given by, 



2mv 0 C 20 
R 2 



Equation 23 
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where M is the magnification of the image associated with the zero order (M 0 ) . Hence, in 
general the magnification of images associated with the non-zero diffracted orders is not equal 
to that of the zero order (for 0 C 20 *0). 

Consider the special case of a single lens and a detector (plane B) placed one focal length from 
the lens so as to form an image of an object at infinity (Figure 19). Adding a quadratically 
distorted grating (4) to the lens (LI) effectively decreases its focal length in the +1 diffraction 
order and increases its focal length in the -1 order. Using standard lens formulae with an image 
distance of one focal length, the object planes imaged in each diffraction order are at locations 
given by, 



Hence, the images corresponding to diffraction orders of equal magnitude but opposite sign (for 
example planes 5 and 7 in Figure 19) are placed symmetrically, equal distances either side of 
the lens, while the zero order is an image of a plane at infinity. The image magnification 
relating to each object plane is, 



The images corresponding to diffraction orders of equal magnitude but opposite sign therefore 
have magnifications of equal magnitude but opposite sign. 

Broadband Operation 

The angle (0 m ) into which light is diffracted from a grating is proportional to the wavelength of 
incident light(equation 1 6). If a broadband illumination source is used the non-zero orders 
become spectrally dispersed (along the x-axis in the image plane) and the images associated 
with these orders become blurred. A wavelength-dependent grating period would offer one 
means of superimposing the images at each wavelength and preventing blur. With a standard 
undistorted grating this is unachievable, but a grating distorted according to a quadratic 
function only offers a route to achieving this. 



u 



m 



2m 0 C 20 



Equation 24 



2mf 0 C 20 
R 2 



Equation 25 
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It was shown in equation 1 8 that a shift in the origin of the quadratic distortion function along 
the x-axis causes the separation of the images in the x-direction to change, whilst leaving the 
defocusing power of the grating unaltered. This property can be thought of as arising from the 
variation in grating period along the x-axis, which is generated when the grating is distorted 
according to a quadratic function. The grating period (d) at a distance x from the origin of such 
a grating is given by, 

d = — Equation 26 

*R 2 -2d 0 0 C 20 x 

where d 0 is the period at the origin of the quadratic distortion function. It is the grating period at 
the centre of the incident beam (the centre of the grating aperture when the full aperture is used) 
which defines the image positions along the x-axis in the non-zero orders. If the grating is 
illuminated with a beam which has a radius less than that of the grating aperture, then as the 
beam centre is moved across the grating to a position x along the x-axis, the separation of the 
images changes in proportion to the inverse of the grating period described by equation 26. 

Consider a defocus grating designed to diffract light of wavelength X, at a specific angle, with 
an undersized beam hitting the centre of the grating, where the grating period is defined as 
d 0 =d x>1 . To equalise the diffraction angle at a second wavelength X2 we require, 

JL 

An =dxi-r- Equation 27 

A,, 

Using equations 26 and 27, the diffraction angles at the two wavelengths can be equalised by 
offsetting the X 2 beam centre a distance x n along the x-axis relative to the X { beam centre 
where, 



x x2 - 0 r - (^2 -M 

Z 0 < -20 Q X1 



Equation 28 
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Note that the required offset is proportional to AX (AX =VA.,). It follows that if the input beam 
is dispersed across the grating at the corrected level of dispersion the images can be 
superimposed at all colours. 

For example, consider a grating of radius (R) 0.5cm, grating period (d xl ) 200jim and oC 2 o = 10^i 
at 550nm. This gives a 1 st order diffraction angle of sin9 = 0.0055 radians when the beam hits 
the grating centrally. From equation 28, to obtain the same diffraction angle at a wavelength of 
500nm requires a beam shift of 0.3 1 3cm, and at 600nm a beam shift of -0.3 1 3cm. Figure 20 
shows the beam positions at 600nm (13), 550nm (14) and 500nm (15) relative to the distorted 
grating (4). The technique can be extended to any number of wavelengths or broadband 
illumination subject to equation 28 being satisfied at all wavelengths. 

In order to implement this technique a dispersive system can be introduced into the optical 
system before the distorted grating. The dispersive system must introduce a beam offset 
proportional to wavelength (equation 28) and leave the propagation direction equal at all 
wavelengths. Figure 21 is a schematic of such a system where broadband radiation 16 which 
enters the dispersive system 17 is linearly translated according to wavelength with 3 discrete 
output wavelengths, 1 8, 19 and 20 from the continuous bandwidth shown. Note that the output 
radiation 18, 19 and 20 all follow parallel paths. 

The construction of a dispersive system with the desired properties is within the knowledge of a 
skilled person and could be implemented, for example, by means of two blazed diffraction 
gratings, projecting the +1 order from the first grating onto the second grating and using the -1 
order from the second grating. Another possible approach is to use a grating/prism or 
prism/prism combination. 

Further information can be obtained from "MODERN OPTICAL ENGINEERING - The 
Design of Optical Systems" by Warren J, Smith, published by McGraw Hill (see especially 
chapter 4.3). 
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Illumination System 

In order to achieve maximum possible resolution and storage density each layer in a multi-layer 
data storage medium must be illuminated with a diffraction limited spot. Figure 22a shows the 
distribution of light that must be achieved from the source 21 . In currently proposed systems 
with a moveable lens (US 5202875), the illumination uses the same optics as the read system. 
The illumination is focused at the right depth but needs additional spherical aberration 
correction. 

The distorted grating described thus far and shown in Figure 1 5, images three on-axis objects 
onto three spatially separated image positions (Figurel4). If the same grating is used in an 
illumination system, with a single illumination source, the spots of illumination are focused on 
the correct layers but are laterally displaced as shown in Figure 22b. Such a grating cannot 
therefore be used for both illumination and reading from a multi-layer structure. 

However, using the techniques described herein, shifting the origin of the quadratic distortion 
function along the x-axis, through the parameter x 0 in equation 15, allows the diffraction orders 
to be aligned along the optical axis. A grating of this type (Figure 1 8), incorporated into optical 
system 1 of Figure 22b would produce the desired illumination shown in Figure 22a. If the 
grating were additionally distorted to include the spherical aberration terms in equation 1 1 then 
the system would illuminate each data layer as shown in Figure 22a with automatic spherical 
aberration correction. 

Complete System 

Using two gratings of the type described in this invention a complete system for illumination 
and reading from a multi-layer optical data storage medium, with no moving parts and 
automatic spherical aberration correction can be described. The system is shown schematically 
in Figure 23. Light from the source 21 passes through a grating 4b of the type described in the 
paragraph above (and Figure 18), which produces multiple on-axis, spherical aberration 
corrected foci corresponding to the different data layers 5, 6 and 7. Light reflected from the data 
layers passes through a grating 4a of the type shown in Figure 1 5c with spherical aberration 
correction, which produces spatially separated images of the different data layers on plane B. It 
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is well known that a polarisation sensitive beamsplitter and polarisation rotating plates can be 
used to minimise losses in such a system. Numeral 1 is used generally to designate an optical 
system. 



In such a system the terms "object plane" and "image plane" can be confusing: layers 5, 6 and 7 
are "image planes" according to the invention with respect to an "object plane" containing 
illumination source 21 and are "object planes" with respect to "image plane" B. For 
convenience we say that the "image planes" of one apparatus according to the invention are 
coincident with the "object planes" of the other. 

Other embodiments of the complete system will be obvious to the skilled person: for example 
the illumination source 21, grating 4b and associated optical system 1 could be located on the 
same axis as the imaging plane B and grating 4a, located on the other side of 'object' planes 5, 
6 and 7. In such an embodiment, planes 5, 6 and 7 are illuminated in transmission. 
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Claims 

1. An apparatus for producing simultaneously a plurality of spatially separated images from an 
object field comprising: 

an optical system arranged to produce an image associated with a first focus condition; 

a diffraction grating arranged to produce, in concert with the optical system, images associated 
with each diffraction order and 

means for detecting the images, 

wherein the optical system, diffraction grating and detecting means are located on an optical 
axis and the diffraction grating is located in a suitable grating plane and is distorted 
substantially according to a quadratic function so as to cause the images to be formed under 
various focus conditions. 

2. The apparatus of claim 1 where the function according to which the grating is distorted 
includes further terms for producing different amounts of spherical aberration in the images 
associated with each diffraction order. 

3. The apparatus of claim 2 where the spherical aberration of images associated with each 
diffraction order is arranged to correct for spherical aberration associated with the different 
depths of substantially parallel planes in object or image space. 

4. An apparatus according to claim 1, 2 or 3 whereby the origin of the distortion function of the 
diffraction grating is displaced from the optical axis. 

5. An apparatus according to claim 4, whereby the origin of the quadratic distortion function is 
displaced to cause alignment along the optical axis of the images associated with each 
diffraction order. 
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6. An apparatus according to claim 1 whereby the diffraction grating comprises a set of two or 
more diffraction gratings designed such that the various diffraction orders are spatially 
separated. 



7. An apparatus according to claim 1 whereby the diffraction grating is any one of an 
amplitude-only diffraction grating, a phase only diffraction grating or a phase and amplitude 
diffraction grating. 

8. An apparatus according to claim 1 whereby the diffraction grating is polarisation sensitive. 

9. An apparatus according to claims 1 whereby the diffraction grating comprises two gratings 
sensitive to different polarisations and arranged such that the diffraction orders produced by 
said gratings are spatially separated. 



10. An apparatus according to claim 1 whereby the diffraction grating is a programmable 
grating. 



1 1. An apparatus according to claim 1 whereby the diffraction grating is a reflective grating or a 
transmissive grating. 



12. An apparatus according to claim 1 whereby the grating is any of a two-level (binary) 
structure, a multi-level (digitised) structure or a continuous-level (analogue) structure. 

13. The apparatus of claim 1 adapted for forming images on a plurality of image planes, from a 
single object plane. 

14. The apparatus of claim 13 where the object plane contains a source of illumination which is 
used to illuminate the image planes. 

15. The apparatus of claim 1 adapted for producing substantially in focus images in a common 
image plane, from a plurality of object planes. 
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16. The apparatus of claim 15 and further including apparatus according to claim 14 wherein 
the object planes of claim 15 are coincident with the image planes of claim 15. 

17. The apparatus of claims 15 or 16 where each object plane contains an array of elements, 
capable of existing in at least two states and in which the detector means is capable of 
distinguishing between said states. 



18. The apparatus of claim 1 7 adapted for reading data from a three dimensional optical storage 
medium wherein object planes are located within the medium and the detecting means is 
capable of producing a signal dependent on the state of the elements. 

19. An apparatus according to claim 1 and further including a dispersive system for introducing 
an offset to an input beam of radiation, said offset being perpendicular to the optical axis and 
proportional to the wavelength of the input radiation, whilst leaving the beams at each 
wavelength following parallel paths. 
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